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Various notions from geomètric control theory are used to characterize the behavior of the Marko- 
vian master equation for N-level quantum mechanical systems driven by unitary control and to 
describe the structure of the sets of reachable states. It is shown that the system can be accessible 
fSJ ■ but neither small-time controllable nor controllable in finite time. In particular, if the generators of 

quantum dynamical semigroups are unital, then the reachable sets admit easy characterizations as 
they monotonically grow in time. The two level case is treated in detail. 
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INTRODUCTION 



The main question that we discuss in this work is the following: to which density operators can we drive the 
quantum Markovian master equation by means of coherent control? This problem is of relevance whenever 
one is interested in quantum state manipulation in presence of nonunitary evolution, for example in the 
context of quantum information processing [D], [Ï7| and of molecular control pO[ . The ultimate goal is 
obviously to know when and how the state of a quantum mechanical system can be arbitrarily manipulated 
by means of unitary (reversible) control operations or at least to what extent this is possible. 

The viewpoint we take in this work is that of "classical" geomètric control theory which provides us the 
tools to mathematically formalize and answer the qüestions posed. In classical control terms, the set of 
density operators to which we can steer the system is called the reachable set and the problem of arbitrary 
Oh, manipulability of the state can be formulated as a controllability problem. 

The infinitesimal structure of the so-called quantum Markovian master equation, i.e. the "axiomàtic" 
model for an open quantum system, is known since the works of Lindblad JÏ4| and Gorini-Kossakowski- 
Sudarshan and it is a prerequisite for the utilization of the Lie algebraic controllability methods developed 
below. We use the so-called vector of coherences formulation for the density matrix EJ, i.e. the expectation 
vàlues corresponding to a complete set of Hermitian operators, here the Gell-Mann matrices. Such formu- 
lation allows to treat the master equation as a control system with affine vector fields or, geometrically, as 
a system living on a homogeneous space of a matrix Lie group and subordinated to an affine group action, 
plus constraints originating from the complete positivity of the quantum dynamical semigroup. If we drop 
these constraints, the system falls into a class of systems whose controllability properties were studied in 
detail in the past, see || [ï^] for a general overview, |], ^, [Ï0| for the particular case of affine fields. Including 
the complete positivity requirements totally alter these results, because of the relaxation it induces. 

The qualitative difference between studying the master equation and its controllcd counterpart is that the 
master equation is an ODE whose solution, obtained integrating a single vector field, is a one parameter 
semigroup; the presence of control parameters in the controlled master equation implies that we have to 
consider a family of vector fields simultaneously, and therefore the admissible flow is a multiparameter 
semigroup or Lie semigroups [||. Such semigroup is the reachable set. When the reachable set is large 
cnough to be a subgroup or at least to act transitively on the homogeneous space, then we have controllability. 
The problems arise when the reachable semigroup is not a group, as in the case of the controlled master 
equation. A novel element with respect to, for example, the control of Schròdinger equations (^] is that in the 
master equation one has to deal with a true drift term, i.e. a vector field which is both noncontrollable and 
nonrecurrent. Then it can happen that although it is (often) possible to generate motion in any direction 
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(i.e. we have the accessibility property), the system in never controllable in finit e time because the flow 
cannot be reversed. In other terms, the reachable set may be opcn and dcnse in the spacc of admissible 
density operators, but the initial condition of the controlled master equation always lies on the boundary of 
such set for any finite time and therefore it is not possible to reach arbitrary points in its neighborhoods. 
The vector of coherences representation is very useful in this respect, as it allows to explain the lack of 
controllability in terms of the trace of the dissipation/relaxation superoperator. In fact, the main reason for 
noncontrollability lies in the structure of the nonunitary operators given by complete positivity. When such 
infinitesimal gencrators is unital this is clcarly visible: for the density operator p, tr (p) gives the level sets 
of a quadratic Lyapunov function centered in the origin. In this case, the controlled dynamics is stable and 
the control alone allows only to move within one of the level sets, not to pass from one level set to another. 
Since the nonunitary operator is pointing inward, as time passes also the controlled integral curves can move 
only inward and this establishes a monotonicity relation among the sets reachable at different time instants. 
As pointed out for example in |2p|] , the presence of a dissipation (nonunitary) operator is essential for any 
motion not confined to a sphere in l n to be accomplished. Notice that this holds regardless of the existence 
of a thermodynamic equilibrium, i.e. a fixed point for the original uncontrolled master equation. For afBne 
dissipation operators, the situation is slightly morc complicatcd and controllability may be recovered as a 
limit process. The atom with spontaneous emission is one such case and will be discussed in some detail. In 
this case, motion is not confined to the inward of spheres in R™ and "purification" processes are possible. 
The organization of the paper is as follows: in Section we review all the relevant notions concerning 



III the formalism of 



controllability of bilinear/affine systems on homogeneous spaces of a Lie group; in Section ' 
the vector of coherences parameterization is recalled and used to discuss controllability of Liouville dynamics; 
in Section MW the controllability of the master equation is treated and the main Theorem formulated. Finally, 
in Section jyTthe two- level case is discussed in detail, first for general dissipation operators and then for few 
significant examples. 

It is worth remarking that all our considerations make sense for finite dimensional quantum systems. 



II. DRIFT AND CONTROLLABILITY FOR BILINEAR CONTROL SYSTEMS 

All properties introduced in this Section are Standard in geomètric control and are adequately surveyed 
for instance in fi |Ï8). Consider the following bilinear control system: 

i 

x = B x + UkB k x 

k=l W 

x(0) = Xi 

where the controls u±, . . . , u q are real valued piecewise constant functions defined on [0, oo), Bq, . . . ,B q are 
square matrices and x £ M, an analytic manifold of real dimension n. In this work: M is Mg = M. n \ {0} or 
some n-dimensional homogeneous space (like a sphere) contained in W , r > n, or some subset of M™ like a 
sòlid unit ball. The vector field Bqx is called the drift and BiU\x, . . . , B q u q x are the control vector fields. 
Given Xi £ M, let us call 1Z(xí, T) the reachable set from xt at time T > for the system (Q) : 

1Z(xí, T) — {x £ M s.t. x(0) — Xi and x(T) = x, T > 0, for some admissible control U\, . . . , u q } 

If 1Z(xí, < T) = Uo<t<T'R-( x ii t), then the reachable set from Xi is IZ(xi) = Uck^og^^íi t). In correspon- 
dence of a given Tt and of TZ(x{, < Ty), one can define the notions of finite time controllability as follows: 

Definition 1 Given Tf > 0, the system (|ï|) is Ty-controllable i/V Xi, xj £ M 3 admissible control functions 
U\, . . .u q such that the flow of (|ï|) satisfies x(0) = Xi and x(Tf) = Xf. 

The existence of a Tf finite is important for the application discussed in this work. When, instead, we are 
interested in controllability for any time in 1Z(xí), then we can use the following: 

Defïnition 2 The system ([!]) is controllable if any Xf £ M is reachable from any Xi £ M for some admissible 
control function u±, . . .U q . 
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Unlike the reachable set which accounts only for the positive time evolution of the trajectòries of the 
system, the orbit 0{xí) requires to consider complete vector fields, i.e. defined on the whole time axis: 

0{xí) = UígK {x £ M s.t. x(0) = Xi and x(t) = x, t £ R, for some admissible control u\, . . . , u q } 

The difference between 0(xí) and 7Z(xí) corresponds to the difference between the accessibility and con- 
trollability properties. 

Definition 3 The system ([ï]) is accessible iflZ(xi, < T) contains nonempty open sets of M for all T > 0. 

While accessibility guarantees the existence of open reachable sets, it does not say anything on ce, belonging 
to it. 

Definition 4 The system (Q) is small-time controllable if Xi belongs to the interior of the reachable set, 
mtH(xi, T), for allT > 0. 

The accessibility property admits an algebraic characterization in terms of the Lie àlgebra generated by 
the vector fields Bqx, B\x, . . . , B q x, call it Lie (Bqx, B\x, . . . B q x). 

Theorem 1 (Lie algebraic rank condition (LARC)) The system is accessible if and only if 

dim (Lie (Bqx, Bix, . . . B q x)) = dim (M). 

For bilinear systems, when accessibility holds there exists a Lie group of transformations, call it G, of 
(finite) dimension greater or equal than M acting transitively on M and to which we can lift the system. 
Invariance of the vector fields on a Lie group implies that the controllability conditions are global and 
independent of the point of application. For example, for both TZ{xí) and 0(xí) we have TZ(xí) — TZqXí and 
Oixi) = OqXí, with TZg — T?-(I) and Og = O(I) reachable set and orbit of the lifted system, where / is the 
identity matrix of G. Therefore we can work indiffcrcntly with vector fields on M (Bqx, B\U\x, . . . B q u q x) 
or with right invariant vector fields on the Lie group G (i.e. the matrices Bq,BiU\, . . . B q u q of g, the Lie 
àlgebra of G), to which we have lifted the system, starting from the identity of G: 



In particular, the LARC condition and the so-called orbit theorem guarantee that Og is the whole G and that 
M is nothing but a homogeneous space of G expressed in terms of equivalence classes as Gx, x £ M. The 
Lie àlgebra g is therefore equal to Lie (B , . . . B q ) and the accessibility condition reformulates as transitivity 
of G (or of g, with a common abuse of terminology) on M. 

Theorem 2 The system (^) is accessible if and only if g is transitive on M. 

The LARC condition is only a necessary condition for controllability, even when it holds the reachable set 
needs not be the whole Lie group G. When TZg Ç Og the lifted system is not controllable, the reason 
being that the drift is allowed to flow only along the time forward direction and may not be reversible by 
means of the control vector fields. In fact, the control vector fields are "complete" in the sense that, since Uk 
can assume both positive and negative vàlues, once exponentiated they generate a one parameter subgroup 
exp (tukBk, t>0). On the contrary, the drift produces only a subsemigroup exp (tB , t > 0) and thus TZg 
in general only has the structure of a Lie semigroup of G H. 

The case of a compact group is exceptional, since compact Lie groups do not admit semigroups: 
exp (tB, t > 0) = exp (tB, t £ R). Hence TZg collapses in Oq and the accessibility property collapses into 
( "long time" ) controllability. The LARC condition then becomes necessary and sufficient for controllability: 
TZ(xí) = O g xí = Gxí = M, Wx t £ M. 

In general, however, one has to deal with the case of TZg being only a Lie semigroup. Even if TZg is a 
proper semigroup, TZg Ç G, it may still happen that the action of TZg on M is transitive. In the literature, 
most results are in the form of sufficient conditions for controllability. For the system (]]]), examples are 

1. TZg — G and G acts transitively on M; 

2. TZq ac ts transitively on M; 

3. xí £ intn(xi, T) VT > 

For our case, none of these (or similar) conditions hold and "negative" results have to be established. 
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A. Afflne vector flelds case 

The case of affine vector fields generalizes (|j) to the following set of ODEs: 



x = B x + b Q x Q + / , {ukB k x + ò fc Xo) 
x(0) = Xi 



9 

k=l 



(3) 



where xq is a real constant. It corresponds to a Lie group of transformations having the structure of a 
semidirect product K(§)V with V typically a n-dimensional real vector space and K a Lie group acting 
linearly on it. The dimension of K(§)V is dim (K ) + n. By choosing the following homogeneous coordinates 
for the state x — [xq, x T ] T , the system (||) recovers the linear form of ([[]): 



1 

x = Bqx + UkBkX 



where Bh 





bk Bk 



. The homogeneous coordinates allow to transform the affine action of K(g)V on x into 
linear action on x. If g = ^ ? G G = T \ 2- 1 the action $ : G x M — > M is: 



£0 



"1 0" 


e G = 


"1 


0" 


v k 


V 


A' 



$(fli)(x) =gx = 
so that the affine vector fields induced on i by $ are: 

$*(B)(x) = Bï = 

and the Lie bracket is: 

[I, S] 





Bx + bxo 





Ab - Ba [A, B] 



Also for affine systems, special sufficient conditions for accessibility and controllability have been devised, 
see §, g, for details. 

III. CONTROLLABILITY OF HAMILTONIAN DYNAMICS 

To describe differential equations for density operators, we make use of the so-called vector of coherences 
formulation. A few essential facts about it are reported below; see for examplc [0 for a thorough description 
and further references. 



A. Density operators and vectors of coherences 

The state of a quantum mechanical system in an A-dimensional Hilbert space Ti. N can be described in 
terms of a trace 1 positive semidefinite Hermitian operator p called the density matrix. If the density 
operator is entirely characterized by a wavefunction \ip), then the system is said to be in a pure state, p is 
dcfined as p{t) — IV'XV'I and tr(p 2 ) = 1. If instead we have a statistical ensemble p(t) — J^iLiP^l^^i^^l 
for pW > and ^.jp!*) = 1, then the system is said to be in a mixed state characterized by the pairs 
{pW, and tr(p 2 ) < 1. In both cases, the properties of Hermitianity p = p^ and of unit trace tr(p) = 1 

imply that the N x N matrix representing the density operator depends on n = N 2 — 1 real parameters. Up 
to the imaginary unit, N x N traceless Hermitian matrices form the Lie àlgebra su(A) of dimension exactly 
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n. If to it we add the (properly normalized) unit vector, then we obtain a complete basis for the density 
operator of an iV-dimensional quantum mechanical system. In fact, the iV-dimensional Pauli matrices Xj, 



see for example [13 for their explicit expression, and the identity matrix Xq = N ^I, form a complete 
orthonormal set of basis operators for p (orthonormal in the sense that tr(AjAfc) = 5jk)- In particular, 
then, p = Y^j=o tr(pAj)Aj = Y^j=o Pj\i w ith p = fixed constant and the n real parameters pj giving 

the parameterization of p. Since the X k , k — l,...,n, form a compete set of observable operators, the 
Pj = tr(pXj) are expectation vàlues of p. Call p = [pi-. ■ p n ) T such vector of coherences of p. Due to the 
constant component along Ao, p is living on an affine space characterized by the extra fixed coordinatc 
po = N~ï . Such n dimensional Liouville space of vectors p — [po p% ■ ■ ■ p n ] T — [pa P T ] T has Euclidean inner 
product given by the trace mètric: ||p|| = \J ((p, pj) = ^/tr (p 2 )- The condition tr(p 2 ) < 1 then translates in 
p-space as p belonging to the sòlid ball of radius 1 centered at [po . . . 0] T , call it B ra , for all positive times. 



The surface of such ball generalizes the idea of Block sphere to N dimension and corresponds to pure states 
||p|| 2 = 1. In tcrms of vector of coherences, the condition tr (p 2 ) < 1 becomes the smaller ball ||p|| 2 < 1— j? 
(centered at the origin). 



B. Hamiltonian dynamics 

If H is a constant fmite dimensional Hamiltonian, for the density operator the Liouville equation is given 

by 

p(t) = -i[H, p] = -i&d H (p) 

If — iH € su(iV), then — iad^f is a so-called commutator superoperator i.e. a linear operator in the n 
dimensional Liouville space of p vectors. In terms of p, the action of — iad# is linear: 

p = -ia,d H p (4) 

H being traceless and Hermitian, in the {Xj} basis: H = hiXi. The process of passing from p to p 

is mathematically equivalent to passing to the adjoint representation of the Lie àlgebra su(N). In fact, the 
corresponding basis in the adjoint representation is given by the n x n matrices ad>j , . . . ad^ of elements 
(adA;)^. = ifijk with fijk real fully skew-symmetric (with respect to the permutation of any pair of indexes) 
tensor. Thus — iadn — —i SILi h>l a >d\ l . The n x n matrices — iad^ , . . . , — iadA„ are real and skew-symmetric 
and as such they are part of a basis of so(n). Since dim(so(n)) = "^" 2 ~ 1 · ) = N =4p +2 , for N > 2 the 
n matrices —iad^, . . . , —iad\ n span only a proper subalgebra ad su (jvi of so(n). For example for N = 3, 
n = dim (su(3)) = 8, while dim (so(8)) = 28! In the Liouville equation (||), the propagator for p corresponding 
to the Hamiltonian H is an orthogonal matrix: 

p{t) = g(t)p(0), g(t) G exp(ad fluW ) C SO(n) 

such that g(t) = —íadngit), g(0) — I. The action of any g € SO(n) on p is isometric and as such it preserves 
the inner product ||p||. 



C. Coherent control of Hamiltonian dynamics 

Assume that the Hamiltonian H is composed of a time-invariant part Hq representing the free evolution 
of the system plus q time-varying forcing terms representing the interaction with q external fields, modeled 
semiclassically: 

q 

H(t) = H Q + J2 u k {t)H k , -íHq, -iH k e su(N) 

with the parameters ut representing the control fields. Consider a pure state of ket \ip) G § C TÍ N (where S 
is the sphere in TV dimensional Hilbert space) and its Schròdinger equation 

i^) = H \^) +J2u k H k \^) MO)) = |Vo) (5) 
fe=i 



G 



The spherc § in Ti N is a homogeneous space of SU(N). Compactness of SU(N) plus transitivity of the 
SU(N) action on S in this case guarantee the following (see jn] for thc original formulation or H for a 
thorough discussion) : 

Theorem 3 The system (^) is controllable if and only if Lie (—íHq, ~iU\, . . . , —iH q ) = su(N) 

By computing the (real) dimension of such Lie àlgebra: dim (Lie (—iHo, . . . , —iH q )) — dim (su(iV)) = n = 
dim(S). 

In the following we will always consider the controllable case for the wavefunction 
Assumption Al The system (|5|) is controllable. 

Passing to density matrices, for a mixed state p driven by the same Hamiltonian H(t) the corresponding 
forced Liouville equation written in terms of vector of coherences is: 

p = -i&d Ho p - i 22 u k &d Hk p (6) 

k=l 

The vector fields — iadjj D , . . . — iadn corresponding to thc Hamiltonian dynamics lack the translation compo- 
nent and belong to a subalgebra ad 5u (jv) of so(n). Just like the Lie group SU(N) is acting transitively on the 
unit sphere on 7í N , so the orthogonal group SO(n) is acting transitively on a sphere ||p|| 2 = const < 1 — ■L·. 
By dimension counting, exp(adwjv)) is not acting transitively on such sphere if N > 2. In fact, it is wcll- 
known that coherent control cannot modify the eigenvalues of p, and so controllability can occur only insidc 
the leaf of the foliation of p (determined by the eigenvalues) that one starts with. See |Ï9| for a description 
of the kinematic equivalence classes of density matrices in the context of dynamical control, or |ï^ | for a 
complete description of the invariants of motion. 



IV. CONTROLLABILITY OF MARKOVIAN MASTER EQUATIONS 

The requirement of tr (p 2 ) < 1 for the density operator is reformulated in the vector of coherences param- 
eterization as ||p|| 2 < 1. Thus B™ has to be made invariant by the quantum dynamical evolution. The main 
feature of the master equation is to capture all the possible infinitesimal generators that fulfill this condition. 
Obviously, also the driven master equation has to live on B n , and all controllability qüestions have to be 
restricted to B". 



A. Master equation 



Calling Cd the superoperator representing the relaxing/dissipating part of the dynamics, in the basis {A^} 
of traceless Hermitian matrices the Markovian master equation is expressed as Q: 

1 - 

p = C H {p) + C D (p) = -iad H (p) + - ^ a jk ([Xj, pX k ] + [Xjp, X k \) 

(7) 



1 

= -i{H, a jk (2X 3 pX k - {X k Xj, p}) 



2 

where the Hermitian matrix A — (ajk) is positive semidefinite, A > 0, and { • , ■ } is the anticommutator. 
For the basis {Aj}, unlike a Lie bracket which is linear in the generators, the anticommutator has an affine 

structure: {Aj, Afc} = ^j^-SjkXç, + X)"=i djklXi, with djki the real and fully symmetric tensor. The expression 
of (Q) in tcrms of vector of coherences is as follows: 

n 

p = -iaà H p + ^ a Jk (LjkP + Vjkpo) (8) 
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with Ljk n x n complex matrix of mixed symmetry and Vj k imaginary n vector given by: 

1 ™ 

-Ljk — i^jk^lr — T ^ > ((/jmr H~~ ^jmr^jfkrnl [fkmr ^jfemf) fjml ) 



(9) 



B. Coherent control of master equations 

Under the assumption of weak and high frequency control fields, it is acceptable to assume that no time 
dependence is induced in the Cd term by the external fields. Adding the controls, eq. (||) modifies as [EM : 



p = C Ho p + ^ Uk.CH k P + C D p 



k=l 



■i&à Ho p - i ^ u k ad Hk p + a,j k (LjkP + VjkPo) 



fc=l 



(10) 



or, in homogcncous coordinates and calling Lj k = 







, j,k = 1, 



p = Cp = C Ho p + 2J u k £ Hk p + £ D p 
k=l 

i 





— iadíf 



fc=i 





— iadjj. 



pel" 

n 

P+ X! ajkLjkP 

3,k=l 



(11) 



The state of ( |ïï| ) is living on R n+1 and is constrained by the positivity of p requirement to belong to B™. 
However, the dissipation term Cd is not coherent and as such it enlarges the integral group of ( |l l·| ) from 
exp (ad su (jv)) to one of the Lie groups properly containing it. Examples are SO(n, R), SL(n, R), GL~(n, K), 
the connected component of GL(n, R) containing the identity, or their semidirect extensions SO(n, R)(S)R™, 
SL(n, R)(s)R n and GL + (n, R)®R n . 

Since A is Hermitian, the number of independent parameters cijk is n 2 . It is convenient to rear- 
range the n 2 degrees of freedom in the following manner. From (^), it is straightforward to check that 



Re [(Lkj)ir] = Re [(Lj k )i r ], Im [(L k j)i r ] = -Im [(L ]k )ir] and therefore L kj = L 



If we call L% 



Re[(L jk ) 



and Lf k — Im [(Lj k )i r ], then with respect to index permutation L^ k is symmetric, L ki = L^ fc while Lj k is 



skew-symmetric, L 



and a 



Im[a 



kj 

jk\l a fcj 



Similarly, the Hcrmitianity of A implies a k j 



djkLjk + a k jLkj — (2 — 5 J 7 £ )a* í fe 





^ 



+ 2a 







J kj 

a* k or, if we write af k 



l<j<k<n 



Reia 



(12) 



> 



To have A > 0, a number of constraints among the aj k must be imposed like, for example, ajj — 
and üjjükk > (af k ) 2 + (af k ) 2 or \a jk \ < (ajj + a kk )/2. 

Our aim here is to draw conclusions about which p can be reached by means of coherent control. In (|ï^) , 
unlike Ylk=i u kCH k , both Ch and Cd have integral curves that can flow only along the positive semiorbit 
and, in control terms, Ch plays the role of the drift and Cd that of a disturbance. Classically, a disturbance 
can be treated for example as a parametric family of vector fields with parameters belonging to admissible 
intervals (here aj k such that A > 0). However, in the case of Cd parametric the master equation becomes 
a differential inclusion and little can be said about its controllability properties. Therefore in this work we 
will assume to be dealing only with a precisely known value of A of hence of Cd ■ 



Assumption A2 The parameters ajk, j, k — 1, . . . , n are fixed and known exactly. 
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Consequently we can treat Cd as a part of the drift term (together with Ch ) and use the tools of Section ||. 

Under the assumption of unitary controllability, the Lie àlgebra g of interest here is the smallest Lie àlgebra 
of real matrices containing ad su (jy) an d Ch + Cd and closed with respect to matrix commutation: 

= Lie {CH a + Cd, Ch 1 , ■ ■ ■ , C~H q ) 
Once q is computed, the system can be lifted to G: 

q 

g = C Ho g + ^2u k C Hk g + £ D g (13) 
fc=i 

The following Theorem gathers various results about accessibility and controllability for system jïï|). 
Concerning controllability, while for Cd unital the results are sharp (and negative), the case of Cd affine 
is more difficult to treat. In fact, in this case, in spite of the lack of small-time controllability it may 
happen that points that are not reachable in short time are reachable for T large enough and even that 
cl (1Z(xí)) = B n asy mpto tically (cl( • ) means closure). The atom with spontaneous emission discussed in the 



examples of Section VA is one such case. Essentially this fact depends on the existence of a fixed point for 
the master equation and on it being on the boundary of B™, <9B n . However, even in this case dM n is reached 
only asymptotically and therefore the system fails to be controllable in finitc time. 

Theorem 4 Under Assumptions Al and A2, we have the following: 

1. the system ( p"l·| ) is accessible in B" if and only if Q = g[(n,R) or q = g[(n,R)(S)R™ 

2. the system ( p"l·| ) is never small-time controllable in B" for Cd ^ 0. 

3. the system ( pi] ) is never Tf -controllable in B" for any Tf > and Cd =/= 0. 

4. the system ( p"l| ) is never controllable in B™ for Cd 7^ unital. 
Proof. 

The proof of Part 1 follows from Theorem ||. The only transitive Lie algebras on Rq (and thus on B n ) are 
sl(n, R) and gí(n, R) and their semidirect extensions sl(n, R)(S)R™ and $l(n, R)(S)R™. Recali that matrices 
in sí(n, R) are traceless and that, using the decomposition gl(n, R) = sí(n, R) © span(/„) (/„ the n x n 

identity matrix), if íï{Cd) = na and ï — ? P , Cd can be split as Cd — cel + Cd, cu G R, a < 0, 

Cd G sí(n, R)(S)R™. But since üjj > 0, in order for Cd to be traceless it must be üjj = Vj = l,...,n and 
hence, from |ajfc| < (üjj +afcfc)/2 all ajk — 0. Therefore only QÍ(n, R) and gl(n, R)(S)R™ are compatible with 
,1 • 0. 

To prové Part 2, one needs to show that the initial condition pi does not lie in int (lZ(p~i)). It is quite easy to 
vcrify it for Cd unital. In fact, if the initial state pi is such that < \\pi\\ = 5 < 1, then for the Hamiltonian 
part (((Cho + J2t=i u kC~H k ) Pi, Pi}} = while Cd is pointing inward: ({CdPí, Pí}} < 0. Therefore, the ball 
of radius S is invariant for the flow of ( |TÏ| ) and pi lies on the boundary of TZ(pi). For Cd affine, the lack 
of small-time local controllability is automatically verificd for pure states \\pi\\ = 1, because the physics of 
the problem imposes that p such that ||p| = \J\ + e, e > is not admissible. Writing the integral curves of 

the control system as p(t) — $ ^Texp / *(£(r)(ir^ (p~i) — g{t)p~i with T the Dyson operator, we can lift the 

dynamics to the system ( |Ï3| ) with initial condition g(0) = I. pi int (TZ(pi)) for = 1 implies that the 
reachable set TZ(pi) — TZgPí cannot be transitive on any neighborhood of pi and that for the lifted system 
/ ^ int^c)- But, due to right invariance, the propcrties of accessibility, controllability and transitivity 
have a global character and therefore IZq is not transitive for any neighborhood of any pi € B™. 

Concerning Part 3, if a finite time Tf is fixed, the reachable set TZ(pi, < Tf) for the master equation is 
always only a Lie semigroup. In fact, if the fixed point of Cd (when it exists) belongs to int (lZ(p~i, < Tf)) then 
c\(K(pi, < T f )) Ç c\(Jl(pi)) Ç B"; if instead it belongs to <9B" then c\(Jl(pi, < Tj)) C cl(^(pi)) = B". 
Even if a fixed point does not exist, we have that the norm of pi can grow only if ((CdPí, Pí}} > and that 
|jp(í)|j can approach 1 at most as t — > 00. Since Ch & n d Ch 1 , ■ ■ -Ch preserve the length, excluding the 
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trivial cases the control cannot speed up the convergence to <9B n from its "best" initial condition. But even 
in that case convergence is only asymptotic. Therefore, for any fixed Tf the open set lZ(p~i, < Tj) cannot be 
all of B" and neither can its closure. 

Finally, the proof of noncontrollability in S™ for Cd unital follows from the same argument used above in 
Part 2. 

□ 

For the system lifted to its integral group, the small-time controllability property collapses into controlla- 
bility and we have the following. 

Corollary 1 The "lifted" system @ is accessible for G = GL+{n, R) or G = GL+(n, R)©R n but it is 
never controllable on G for Cd ^ 0. 

Proof. The first part is obvious, since accessibility on the Lie group is a necessary condition for accessibility 
on the homogeneous space. Concerning controllability, from the proof of Theorem ^, Part 2, for the system 
(||), I i int(ft G ). But, for Lie groups, such property is a global one and therefore the systems in never 
controllable. □ 

Another way to prové the previous Corollary is via piecewise constant controls: in this case, 
g(t) = Texp f Q C(r)dT = H j=1 exp ((C Ha + C D + YX=i u k 3 X-B k ) (tj ~ tj-i)) and using the for- 
mula det(exp(-)) = exp(tr(-)) we have det (jT j=1 exp ((C Ho + C D + J2l=i u kjC Hk ) (tj - = 

ex P (Sj=i tr ((£h + Cd + J2l=i u k 3 £l·l k ) (tj - íj-i))) = exp(tr(£ D )í) < 1. Therefore, g(t) G 

GL + (n, R) or g(t) G GL + (n, R)(S)R n (only in these cases (|ï|) is accessible) must be such that àet(g(i)) < 1 
and cannot generate the whole Lie group. 

Yet another method to show the same thing is to use the necessary condition of Lemma 6.8 of [0. Write 
C D = cd + C D , a e R, a < 0, C D G sl(n, R)©R™. Since [I, F] = VF G gí(n, R)(S)R", under the 
assumption of accessibility the ideal in g generated by the control vector fields coincides with the derived 
subalgebra Qi = [g, g) (and span(J) with the center of g) and it is contained in si(n, R)(S)R Ii . A necessary 
condition for TZq = G is that exp ((Ch + E-d) t) G exp (gi) (i.e. to SL(n, R) or to SL(n, R)(S)R") for some 
t > which is obviously never true. 

For Cd unital, the reachable sets are balls in R™ (centered in 0) and are completely characterized by the 
following monotonicity property: 

Corollary 2 If the system ( |ï"ï| ) is accessible and if Cd unital then lZ(f>i, < Ti) C lZ(p~i, < T2) V0 < T\ < T2 
and lZ(p~i) is the ball of radius \\pi\\- 

Proof. It follows from the observation above that p u of norm Hp^H lies on the boundary of the set reachable 
from p v by the integral curves of (11). If p v = p(Ti) = $ ^Texp J Q Tl £(r)drj (pi), then lZ(p~i, < T 2 ) 



7Z(Pii < T'í) U lZ(p(Ti), < T2 — Tx). Notice that this does not require Cd to have a fixed point. Accessibility 
of (fi"ï|), in fact, guarantees that p(t) can be placed on any point of the sphere of radius ||p(í)|| and therefore, 
as t ~ + 00, (|rï|) can be made to tend to the origin regardless of the existence of a fixed point for Cd- Hence 
lZ(p~i) is any thing inside the ball of radius □ 



TWO LEVEL SYSTEMS 



For two level systems, p is the usual Bloch vector. Call Xk = -^°~k, k G {0, x, y, z} the rescaled (identity 
and) Pauli matrices. Then in the {Xk} basis 



Poo P01 
Pw Pu 



P0X0 + p x X x + PyXy + p z X z 



V2 



^75 + Pz Px~ ÍPy 
Px +ÍPy - Pz 



(14) 



and p = [p x p y p z ] T ', where p k = tr (pX k ) i.e. p = p x = V% Re[poi], Py = Im[poi] and 

Pz = 75 ^ 00 _ Pl1 ^- In our case ' ^ Ao ' Afe J' = {^i> = V^^-feAo, V j, k G T = {x, y, z}. Similarly to 
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(O), thc Hamiltonian H can be written as 



H = J2V2h k \k = 



kei 



ïh/f, 



K + ih. y -h z 



and, in the adjoint representation, from — ia,dn = (— *ad_f/) = (J2iex hifipm.) 








—h z 


hy 




"o 





" 




' 


f 




"o 


-1 o" 


iadH = 


h z 





-h x 


= h x 








-1 


+ hy 








+ h z 


1 







-hy 


K 










1 







-1 














(15) 



If {x, y, z} — {1, 2, 3}, specifying the coherent controls: h k = (IiQ k + u k ), k = 1, 2, 3, where ho k are the basis 
components of the time-independent free Hamiltonian Hq and u k — u k (t), k = 1, 2, 3, the control parameters 
(some of the ho k or u k may be 0). In the homogeneous coordinates, the vector field for the Hamiltonian 
acquires only a zero translation, and, from (15), the infinitesimal generators of the coherent rotations are: 



Mi = 



"0 








o - 




"0 








0" 




"0 








0" 

























1 










-1 














-1 


, M 2 = 














, M 3 = 





1 








.0 





1 


. 




.0 


-1 





0. 




.0 








0. 



Hence Ch = J2k=i hç> h M k and Cn k = M k , k = 1, 2, 3. In this case ( |Ï2| ) simplifies to: 





2iaf k Vj k (2 - 5 jk )af k L jk 



(16) 



The 9 degrees of freedom of A (constrained by the positive semidcfiniteness requirement) are captured by 
the 9 real parameters (reindexed cardinally) 



{a 4 , a 5 



■ ai2 



i = ía 3 * a 9 a R a 9 a K a K 1 

/ {"'xyi "'xyi "izi u izi "j/z i "j/zï "sa;! "j/y zz j 



In terms of 04, . . . 0,12, the matrix A is: 



a%o 0,4 + «0,5 a@ + ia'j 
a 4 — ia 5 au as + iag 

06 — ÍO-7 0-8 ~ 0-12 



In order to impose the positive semidefiniteness of A, a sufficient condition is that all the principal minors 
have nonnegative determinant, i.c. 



aio > 0, au > 0, 



ai 2 > 



a w au >a\ + a\, 



a 10 ai2 > ai 



anai2 >a\ + a\ 



(17) 
(18) 



aioanai2 - a w (a\ + ag) - ft U (og + af) - ai 2 (a| + a\) + 2a4(o6a 8 + a 7 a 9 ) - 2a 5 (a 6 ag - a 7 as) > 

The infinitesimal generators corresponding to this parameterization are linear combinations of the Lj k . 
Numbering in the same fashion as the a^ k , a^ k parameters, we obtain the 9 linearly independent generators: 



Ma 



'0 0- 

10 

10 





M, 



i (L, 



0" 





-2 



Mb = L 3 



oooo- 

1 


.0100. 



M 7 = i (L xz - L zx ) = 



"0 0" 



2 

.0 0. 
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M 8 = L vz + L z 



oooo- 


1 
10 



L Z y) — 



- 








0" 


-2 























. 












Mi 



' 

0-10 
-1 



, M n 



0' 
0-100 



-1 



M 12 



"0 








0" 





-1 














-1 





.0 












The above expression of the matrix generators is very convenient for our purposes, because it splits 
the affine and linear parts of the action on p. Furthermore, it makes it straightforward to check that 
Lie (Mi, ... , M i2 ) = flí(3, R)(S» 3 (recali that dim (g[(3, R)(S)R 3 ) = 12). 

In terms of the coherence vector and using homogeneous coordinates, eq. ( |ïï| ) becomes: 

3 3 12 

P = ^h Qk M k p + ^u k M k p + ^a k M k p, pel 3 (19) 

fc=l fc=l fc=4 

Given C d , the corresponding Lie àlgebra is 

Í3 12 3 \ 

h 0k M k + aA, Y u k M k Ç fl [(3, R)©R 3 ; 
fe=l k=4 k=l / 

In general, q varies with the vàlues of a k . A few prototypes of subalgebras obtained disregarding the 
assumption A > are reported in Table |. Once A > is imposed, case 1 , case 3 and case 4 are not anymore 
admissible (the argument is the same as in the proof of Theorem ^, Part 1). 



TABLE I: Examples of subalgebras of g[(3, R)(S)R 3 obtained for different A (not necessarily A > 0) 





CoefHcients aj k 





case 1 


a 5 = a 7 = ag = 0; aio, au, ai2 s.t. tr (Y^,k=io a kM k ) = 


s((3, R) 


case 2 


as = a 7 = a 9 = 0; aio, au, ai2 s.t. tr (J2l=io a kMk) 


fl [(3, R) 


case 3 


ai = ae = ag = aio = au = ai2 = 


ad su(2) ©R 3 


case 4 


aio, au, ai2 s.t. tr (J^lLw a kM k ) = 


«1(3, R)©R 3 


case 5 


aio, au, ai2 s.t. tr (J2k=io a kMk) 7^ 


£|[(3, R)(|)R 3 


case 6 


a4 = • • • = ag = 0, aio = au = ai2 


ad 6u(2 ) © span(/) 


case 7 


04 = a6 = as = 0, aio = au = ai2 


(ad ou(2 ) © span(J))©R 3 



The two level version of Theorem [| is then: 
Theorem 5 For a two-level master equation, under Assumptions Al and A2 we have: 

1. the system ( |Ï9| ) is accessible in B 3 for Cd 7^ olÏ + X)fc=5 7 9 a kM k , a G R, a < 0. 

2. the system ( |Ï9| ) is never small-time nor finite-time controllable in B 3 for Ld 7^ 0; 

Proof. The assumption Cd ^ aï + J2 k =5 7 9 a kM k , a < 0, rules out case 6 and case 7 of Table [| By 
exclusion, or by exhaustive computation using the structure constant of Appendix any nonnull Cd such 
that £p ^ aï + J2 k =5 7 9 a kM k generates the Lie àlgebra of case 2 or case 5 as required by Theorem ||. □ 
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A. Examples 

In quantum information processing, some of the admit well-known physical interpretations in terms of 
nonunitary quantum operations on a qubit normally used in the theory of error correction. For example Mio, 
Mn and M12 are, respectively, the infinitesimal generators of the one-parameter semigroups corresponding 
to bit flip, bit-phase flip and phase flip channels (see |ï(|, Sec. 8.4) and so a depolarizing channel has aio, 
au and 012 all nonnull and equal. 

1. Depolarizing channel (Cd unital) 

The depolarizing channel is given by Cd — al, a < 0. Since Cd commutes with everything, in this case 
the system is not accessible and furthermore its integral curves are not at all modified by coherent control. 
They will be pointing "isotropically" to the origin in M 3 . 

2. Phase flip (Cd unital) 

The phase flip channel is also known as phase damping or pure coherence channel and it is given by Cd 
aligned with M12. The cffcct of this one-parameter semigroup is to "contract" the Block sphere along the 
X x and X y directions, leaving it untouched along A z . As an example, check the accessibility property in 
correspondence of the following simple master equation: 

p = (u\Mi + «2M2 + U3M3 + C112M12) p 

i.e. controls available along all the three directions and no free Hamiltonian. The Lie àlgebra 
Lie {Mi, M 2 , M 3 , M12} is computed using the structure constants given in Appendix |A| and the Jacobi 
identity to eliminate terms not linear ly independent. 

• first level Lie brackets 

[Mi, M12] = -M 8 , [M 2 , M 12 ] = M 6 

• second level Lie brackets 

[Mi , M 6 ] = -Mi, [Mi , M 8 ] = M 12 -M11, [M 2 , M 6 ] = M w - M u 

Therefore 

= {Mi, M 2 , M 3 , M A , M e , M 8 , M 10 - M 12 , M n - Mu, M12} = flI(3,K) 

and the process is accessible. Notice that Miq — M12 and M12 — Mu are traceless i.e. they belong to st(3, K) 
(unlikc A/12) and therefore that 0i = st(3,E), as expected. 

In this case, as it is easy to check (see also 0, Part 2, Sec. II. 5) Cd is not uniquely relaxing, i.e. a 
fixed point for the uncontrolled system does not exist. Thus the asymptotic value depends from the initial 
condition and lim t ^oo p — [po p z (0)] T . Once the control is added, however, the controlled system can be 
made to converge to any p z in the interval [— p 2 (0), p z (0)]. In particular, if p z = then lZ(p(0)) = {p S 
B" s.t. ||p|| < ||p(0)||}. 

A Lie algebraic method per se is normally not constructive. However, what it telis in this case is that 
full accessibility is achieved only at the second level of brackets. Therefore a series expansion cannot be 
truncated before that, if one wants to assure the generation of group actions in arbitrary directions. 

To understand what is happcning to the integral curves of the system, it is convenient to split Cd into 
part in sï(3) and part in gl(3) \ sl(3) as in the proof of Theorem ||: Cd — al + Cd, a < 0, Cd G sí(3). For 
Cd = a í2 Mi2 

00" 
0-100 
0-10 
2 



C D = a 12 Mi2 = al + C D = — rp-f + 
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If F(t) = C D + ELi(K + u k (t))M k , since [I, F] = VF e 8 [(3, 
as the exponential 



thc flow of the system can be written 



g(t) = Tcxp / (al + F(t))6t = exp(íal) T exp / F(r)dr 
Jo Jo 

and its action on p as p(t) — exp(íaï)Texp J Q i r (r)drp(0). The "isotropic" contraction exp(toï) corresponds 
to the depolarizing channel a/2(Mio + Aín + AÍ12) and cannot be reversed. Furthermore, the complete 
positivity constraint A > imposes that the al part must be dominant with respect to the F(t) part. Thus, 
regardless of the control action, the overall result is a contraction of the fiow in B 3 . 



3. Amplitude damping (Cd affine) 



In terms of the master equation, the amplitude damping channel corresponds to an atòmic system with 
spontaneous emission. In a two-level system, the excited state |1) can decay to ground state |0) while 
emitting a photon. The process of spontaneous emission is characterized in terms of the atòmic ladder 
operators <r± = a x ± ia y and of the damping coefficient 71 (7^ > 0) as (see ex. the survey [0 and |Ï6|], Sec. 
8.4.1) 



-i&d H (p) 



71 



(2u-pa + — <T + <j-p — po+oJ) 



(20) 



If, for example, H = \[2 í ho 3 A3 + J2k=l u k^k ) then 









































-ho 3 - u 3 











1 

2 








p = 





h 0s + u 3 





-Ui 


P 


-71 








1 

2 










-u 2 


Ul 


_ 






1 








-1 




3 

M 3 p + 53 


u k M k p + 


^ (M10 - 


-M11 - 


M 5 ) 


p 





(21) 



k=l 



Since the unital part of Cd is not proportional to the identity, by Theorem |^ spontaneous emission is an 
accessible process. From (16h, the relaxation matrix is 



A = 2i 



1 -i 
i 1 




As the unital part of Cd is invertible, (|2lJ) has a fixed point. Since A > but not A > 0, Cd lies on an 
exposed face of the cone A > and in correspondence the fixed point lies on cffi 3 , i.e. it is a pure state. 
Thus, asymptotically, Cd admits a reachable set such that cl(7?.(p(0))) = B . As mentioned in the proof 
of Theorem ||, coherent control can speed up the "purification" of p only for certain vàlues of the initial 
condition. In fact, from 



íc=i 

only the last term gives a nontrivial contribution and can become positive for example in correspondence of 
p z positive and p x , p v small enough: 



71 



(((M10 + M n - M 5 ) P, P)) = ^ (pop, - ^ {pI + P \ + P^)) = ^ íp p, - i||p|| 2 ) 



However, the purification process remains an asymptotic process since as soon as coherent control has brought 
p to [0 + \p z \], then purification can occur only because of Cd- Once again, notice how the role of Cd is 
essential in move around in the reachable set. 
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VI. CONCLUSION 



The aim of this work is to shed light on the possibility and limits of coherent control of Markovian 
master equations using Standard tools from geomètric control theory. It turns out that there is a remarkable 
complementarity between admissible quantum dynamical semigroups and controllability: in the vector of 
coherences representation an admissible Cd has to have a nonull component along the nonzero-trace one- 
dimensional vector subspace of the Lie àlgebra ofnxn matrices (or its semidirect extensions) . A component in 
that direction guarantees noncontrollability in small and finite time. In the simple cases of unital dissipation 
operators, the fact that the "uncontrollable" direction has dimension 1 allows to obtain an order relation 
among the sets reachable at different times by means of arbitrary coherent controls. 
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APPENDIX A: STRUCTURE CONSTANTS OF l(3,R)©R 3 



For the real Lie àlgebra gt(3, R) in the basis Mi, . . . , M\i, the structure constants, call them Cj k , are real 
but not totally skew-symmetric: 



„3 _ 6 _ 
c l,2 — c l,4 — 




- r 12 - r 8 - r 1 - r 10 - r 4 
" c l,8 — L l,ll _ c 2,3 — c 2,6 — c 2,8 


— c 2,12 — L 3 : 4 


- r 8 

— c 3,6 


- r 9 - r 4 - r 2 - 

— c 3,7 — c 3,10 — C 4,S " 


1 


r 3 




— r 5 — r 5 — r 5 — r 2 — 

— c 5,10 — c 5,ll — c 6,9 — c 6,12 — 


„5 _ 7 _ 
c 7,8 — c 7 : 10 — 


c 7,12 ~ 


r 1 - r 9 - r 9 - 
c 8,ll — C 9,H — 9,12 — 


1 




C?,3 = 


r 4 - r 5 - r 11 - r 8 - r 8 
c l,6 — c l,7 — G l,8 — G l,12 — C 2A 


_ r 9 _ „11 _ 


" c 2,9 - 


- r 6 - r 10 - r 6 - - 
■ c 2,10 — c 3,4 — c 3,8 — 


1 






7 4 19 
c 3,9 — c 3,ll — c 4,6 — c 4,7 — 


r 7 - r 3 - 
c 4,9 — c 4,ll — 


r 9 - 


r 3 - r 2 - r 1 - - 
c 6,8 — C 6.1Q — c 8,12 — 


1 
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